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Abstract 

We study cosmological perturbations in two-field inflation, allowing for non- 
standard kinetic terms. We calculate analytically the spectra of curvature and 
isocurvature modes at Hubble crossing, up to first order in the slow-roll parameters. 
We also compute numerically the evolution of the curvature and isocurvature modes 
from well within the Hubble radius until the end of inflation. We show explicitly for 
a few examples, including the recently proposed model of 'roulette' inflation, how 
isocurvature perturbations affect significantly the curvature perturbation between 
Hubble crossing and the end of inflation. 
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1 Introduction 



Inflation provides a simple and elegant scenario for the early Universe (see e.g. [T] for a 
recent textbook presentation). Although single field inflation models are perfectly com- 
patible with the present cosmological data, many early universe models based on high 
energy physics, in particular derived from supergravity or string theory, usually involve 
many scalar fields which can have non-standard kinetic terms. This is why multi-field 
inflationary scenarios, where several scalar fields play a dynamical role during infiation, 
have received some attention in the literature (see e.g. [2]-[T3]). However, except for a 
few specific models, the predictions for the spectra of primordial perturbations are, in 
general, a nontrivial task, in contrast with single-field models. 

The main reason is that the curvature (or adiabatic) perturbation, which is generated 
during infiation and eventually observed, can evolve on super-Hubble scales in multi-field 
infiation whereas it remains frozen in single-field infiation. This is due to the presence of 
additional perturbation modes, often called isocurvature (or entropy) modes, correspond- 
ing to relative perturbations between the various scalar fields, which act as a source term 
in the evolution equation for the curvature perturbation. This property was first pointed 
out in [6| in the context of Brans-Dicke infiation. This phenomenon occurs during infia- 
tion and affects the final curvature perturbation at the end of infiation, independently 
whether isocurvature modes survive or not after infiation. 

The purpose of the present work is to study in detail how the isocurvature perturba- 
tions, present during infiation, affect the curvature perturbations, both at Hubble crossing 
and in the subsequent evolution on super-Hubble scales. Since our intention is to stress 
some qualitative properties specific to multi-field infiation, we have chosen to restrict our 
study to the case of two scalar fields. Moreover, we consider models where a non-standard 
kinetic term is allowed for one of the scalar fields. This includes, in particular, scenarios 
motivated by supergravity and string theory, which have been recently proposed (see, e.g., 

m-m)- 

The production of adiabatic and isocurvature modes for two-field infiation with a 
generic potential, in the slow-roll approximation, was studied in [19] where a decompo- 
sition into adiabatic and isocurvature modes was introduced. Models with non-standard 
kinetic terms for infiatons have been studied in the slow-roll approximation in [9] and [TT] . 
and the adiabatic-isocurvature decomposition technique of [19] was later extended to such 
two-field models in [20] and [21]. Very recently, two-field infiation with standard kinetic 
terms was investigated in [22] at next-to-leading order correction in a slow-roll expansion 
and it was also shown that the adiabatic and isocurvature modes at Hubble crossing are 
correlated at first order in slow-roll parameters. In parallel to these analytical studies, a 
numerical study of the evolution of adiabatic and isocurvature was presented in [23] and 
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In the present work, we extend the previous analyses in the following directions. First, 
we present a detailed analysis of the correlation of adiabatic and isocurvature just after 
Hubble crossing, both analytically and numerically. This correlation was in general sup- 
posed to vanish in most previous works, except in the numerical study |23| and the 
analytical work |22| . both in the context of canonical kinetic terms. Here, we extend 
the analysis with non-standard kinetic terms. We compute analytically the spectra and 
correlation at Hubble crossing in the slow-roll approximation, by taking special care of 
the time-dependence shortly after Hubble crossing. 

Second, we study numerically the whole evolution of adiabatic and isocurvature per- 
turbations from within the Hubble radius until the end of inflation. This allows us to go 
beyond the slow-roll approximation which is needed to derive analytical results. Our nu- 
merical study enables us to see precisely how isocurvature perturbations can be transferred 
into adiabatic perturbations during the inflationary phase depending on the background 
trajectory in held space. We illustrate this behaviour by studying numerically three mod- 
els. The last one is the so-called 'roulette' inflation model, which has been proposed 
recently [T4] . 

The plan of the paper is the following. The next section presents the class of models 
we consider and gives the homogeneous equations of motion as well as the equations 
governing the perturbations. The third section is devoted to the study of the perturbations 
from deep inside the Hubble radius until a few e-folds after Hubble crossing. We then 
discuss, in section 4, analytical methods to determine the evolution of the perturbations 
on super-Hubble scales. Section 5 is devoted to the numerical study of the evolution of the 
perturbations, which is compared with the analytical estimates of the previous sections. 
We finally draw our conclusions in the last section. 



2 The model 



In this paper, we study models with two scalar fields, in which one of the scalar fields has 
a non-standard kinetic term, described by an action of the form 



S 



/\#2 1 



,26(0) 



(1) 



where Mp is the reduced Planck mass, Mp = (SttG)"^/^. This type of action usually 
appears when x corresponds to an axionic component. It is also motivated by generalized 
Einstein theories [6l[II]. When 6(0) = one recovers standard kinetic terms for the two 
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fields. In this section, we give the equations of motion for the homogeneous fields and 
then for the linear perturbations, following the results (and notation) of [20] and [2T] . 
where the same type of models was considered. 



2.1 Homogeneous equations 



Let us start with the homogeneous equations of motion. We assume a spatially fiat FLRW 
(Friedmann-Lemaitre-Robertson- Walker) geometry, with metric 



where t is the cosmic time. One can also define the comoving time t = / dt/a{t). 
The equations of motion for the scale factor and the homogeneous fields read 

4> + 3H<P + V^ = h^e^\^ , 



(2) 



and 



H 



2Ml 



1 . p26 



'9 + e^^x' 



(3) 
(4) 
(5) 

(6) 



where H = d/a and a dot stands for a derivative with respect to the cosmic time t and a 
subscript index or x denotes a derivative with respect to the corresponding field. 



It is also useful to introduce the following slow-roll parameters 



v2 
2b A 



-XX 



2MlH^' 



ViJ 



Vij 
3H^ 



and 



+ e 



XX 



H 



(7) 
(8) 

(9) 
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2.2 Linear perturbations 



We now discuss the linear perturbations of our model (one can find a detailed presentation 
of the theory of cosmological perturbations in e.g. [H, [25l [26j and a pedagogical introduc- 
tion in e.g. [27]). For simplicity, we shall directly work in the longitudinal gauge. In the 
absence of anistropic stress (the off-diagonal spatial components of the stress-energy ten- 
sor), which is the case when matter consists of scalar fields, the metric in the longitudinal 
gauge is of the form 

ds^ = -(l + 2$)dt2 + a2(l-2$)dx2. (10) 
where only scalar perturbations are taken into account. 

We now decompose the scalar fields into their homogeneous (background) parts and 
the perturbations: 



0(t,x) = 0(t) + 50(t,x) and x{t,^) = xit) + Sxit,x). 



(11) 



We shall work with the Fourier components of the perturbations, 50k (^) and Sxkit), 
routinely omitting the subscript k to shorten the expressions. The perturbed Klein- 
Gordon equations read 



+ 3H5(I) + 



+ V^^Sx - 2b^e^'xSx 



and 



40$ - 2V^<^ 



Sx + i3H + 2b^<P)Sx + 



(12) 



6x + 2b^xS(p + 



+ 



e-'" {V^^ - 2b^V^) + 2b^^(j)x 



4x$ - 2e-'^%<^ . 



(13) 



The energy and the momentum constraints, given by Einstein equations, are, respec- 
tively: 



1.2 

3ii'($ + i/$) + + — $ 



2M|, 



+ e'^'x 5x + b^e'''r5<P + V^5<P + V^5x 



-,2fe -,2 



2Ml 



+ e^'x5x 



(14) 
(15) 



It is convenient, instead of using the perturbations Sep and 5x^ defined here in the 
longitudinal gauge, to introduce the so-called gauge-invariant Mukhanov-Sasaki variables: 



Q^ = 5(P+ j^^ and Q^ = 5x + j^^, 



(16) 



which can be identified with the scalar field perturbations in the flat gauge. 

Substituting (fT6l ) into (fT2ll -(fT3ll and using the background equations of motion as well 
as the energy and momentum constraints, one finds 

+ ?,HQ^ - 2e^%^ xQx+{^ + C^,^ + C^^Q^ = (17) 

+ 3HQ^ + 2b^ + 2b^ xQ<p+ (^^ + C^x^ + C^^Q^ = , (18) 

where we have defined the following background-dependent coefficients: 

C,, - -2e b,x +M2-^M^2-^M^2-^ + Mp^ + ^^^^ 

^ 3e^Vx e^^0x^ e^^0^x , Wx , ^^'^ , ..^^ 



-2b, 



^^^^MlH^' ^^^^ 



The two equations fflTl) and flTSl) form a closed system for the two gauge-invariant quan- 
tities and Q-^. 



2.3 Decomposition into adiabatic and entropy components 

As originally proposed in [19], in order to facilitate the interpretation of the evolution of 
cosmological perturbations, it can be useful to decompose the scalar field perturbations 
along the two directions respectively parallel and orthogonal to the homogeneous trajec- 
tory in field space. The projection parallel to the trajectory is usually called the adiabatic, 
or curvature, component while the orthogonal projection corresponds to the entropy, or 
isocurvature, component. Note that there was a semantic shift in the terminology since 
one used to call adiabatic and entropy modes during infiation the two particular solutions 
for the perturbations that would match after infiation, respectively, to the adiabatic and 
isocurvature modes defined in the radiation era. This terminology is used for example in 
the papers on double infiation such that [5] and [TO] . 
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This decomposition into instantaneous adiabatic and entropy components, introduced 
in [19], has recently been extended [28] to fully nonlinear perturbations in the context 
of the covariant nonlinear formalism introduced in [29l [30] . Here, we consider only the 
decomposition at the linear level, but since we allow for non-standard kinetic terms, we 
will need to generalize the equations to such a case, as was done in [20]. Let us recall here 
the main results. 

The essential idea is to introduce the linear combinations 

Sa = cos 9 5(f) + sin 9 5x and 5s = — sm9 dcj) + cos9 e'' 5x , (23) 

where 

cos^ = -, sin^ = ^ with &= J^^Ve^. (24) 

The notations & and 6a are just used for convenience; they do not refer to any scalar field 
a. 

Instead of 6a, it is in fact more convenient to work directly with the Mukhanov-Sasaki 
variables and therefore to define 

Qa- = cos9 + sm9 and 6s = — sm9 + cos9 e'' , (25) 

by noting that 

Q^ = 6a+^<!>. (26) 

In the so-called comoving gauge, the perturbation Q^^ is directly related to the three- 
dimensional curvature of the constant time space-like slices. This gives the gauge-invariant 
quantity refered to as the comoving curvature perturbation: 

n = -Q,. (27) 
a 

The perturbation 6s, called the isocurvature perturbation, is automatically gauge-invariant. 
It is sometimes convenient, by analogy with the curvature perturbation, to introduce a 
renormalized entropy perturbation which is defined as 

S = —6s. (28) 
a 

In field space, Q^^ corresponds to perturbations parallel to the velocity vector (0, e^'x), 
while 6s to the orthogonal ones. 

Introducing the adiabatic and entropy "vectors" in field space, respectively 

Ei = {cos9,e-^sm9) , = {- sin 9 , e''' cos 9) , / = {0, x} , (29) 
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one can define various derivatives of the potential with respect to the adiabatic and 
entropy directions. Assuming an implicit summation on the indices / (and J), the first 
order derivatives are defined as 

K = EiVj, V, = ElVj, (30) 

whereas the second order derivatives are 

= EiEiVij. V.s = EiE^Vjj, V,s = ElEiVij. (31) 

By combining the two Klein-Gordon equations for the background fields, ^ and (Jll), 
one gets the background equations of motion along the adiabatic and entropy directions, 
respectively, 

a + ?,Ha + K = 0, (32) 

d = — hA,(y sinO, (33) 
a 

while the equations of motion for the perturbations read: 

+ ?>HQ, + l — + c,AQa + -y5s + a. 5s = (34) 

5s + mss +[ — + c,sUs- — ^g. + CsaQ^ = o , (35) 



with 



C„ = + + 21/ , + + 26<,(s^l/ - c^V.) (37) 

Css = Vss - (^y^ + b^{l + sj)ceV^ + h^clseVs - &\b^^ + hi) (38) 
C„ - + ^ (30) 



a (t2 MlH 



where sq = sin 9 and cq = cos 9. 



The above system consists fl34tl35p of two coupled second order differential equations 
involving only Q^r and 6s. In order to relate these variables to the metric perturbation 
$ defined in the longitudinal gauge, it is useful to use the Poisson-like constraint, which 
follows from the energy and momentum constraints ( fT4ll and (fTSll . 
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where is the comoving energy density and can be expressed as 

em = &Qa+ [3H + ^] aQ, + KQ. + Ws. (41) 



2.4 Perturbation spectra 

The inflationary observables are customarily expressed in terms of power spectra and 
correlation functions. Given their origin as quantum fluctuations, the perturbations can 
be represented as random variables. We introduce the power spectra of the adiabatic and 
entropy perturbations, respectively 

(g.kQ.k') = ^Po.(fc)5(k-k'), (K'^^k') = -^n.(fc)5(k-k'), (42) 
as well as the correlation spectrum 

27r2 

{QA^s^) = -^CQ^ssimk-k'). (43) 



3 Evolution of perturbations inside the Hubble radius 

In order to study the generation of perturbations from vacuum fluctuations, we start, for a 
given comoving wave number k, at an instant (or Xj) during inflation when the physical 
wave number k/a is much bigger than the Hubble parameter H. Our initial conditions 
are given, as usual, by the Minkowski-like vacuum at Xj, 

Qain) ~ = and 5s{Ti) ^ = (44) 

a{Ti)y2k a{Ti)y2k 

for initial adiabatic and isocurvature fluctuations, respectively. These two initial fluctu- 
ations are statistically independent because the corresponding equations of motion are 
decoupled in the limit k ^ aH. 

Although the adiabatic and entropy fluctuations are initially, i.e. deep inside the Hub- 
ble radius, statistically independent, this is, in general, no longer the case at Hubble 
crossing. In the context of two-fleld inflation with canonical kinetic terms, this point has 
been stressed in the numerical analysis of t23j and studied analytically in [22j . 



In the following, we extend the analysis of [22] to non-canonical kinetic terms. 
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3.1 Equations in the slow-roll approximation 



We start with the perturbations and 5s, whose dynamics is described by eqs. ( l34ll and 
(l35ll . Using the conformal time r and introducing the variables 



these equations can be rewritten in the form 



2K 



u„ 



s I 



a 



a 



Us 

Ua + 
Us + 



a 6s, 



,^ 2^ 

—a + a Ca 



a 



Us 



a + a Cs 

a 



u„ 



0, 



(45) 

(46) 
(47) 



where the four coefficients C/j are given in eqs. (I36l) -(l39l) and a prime denotes a derivative 
with respect to the conformal time r. 

Let us now discuss the slow-roll approximation. The only difference with respect to 
the case with canonical kinetic terms will arise from some of the terms depending on the 
derivatives of h. In the slow-roll approximation, one can use the relation 



Vs 



a 



Ht]^s - h^asl 



and the various coefficients in the above system of equations simplify to yield 



dr^ 



2 + 3e 



1 + 2E--5^ + m4 



where the matrices E and M are given by 
E = 
M = 



-77as 

-6e + 3?7ctct Ar], 



r dr 



^sl 
^sl 

S^SgCe 



Ua 



(48) 



(49) 



-A^sl 
-3^ce{l + si) 



with 



(50) 
(51) 

(52) 



In eqs. (1501 1 and (1511 ). we have kept only the terms linear in b^, i.e. proportional to ^, 
and neglected the terms quadratic in as well as the terms proportional to fo^^. We thus 
treat ^ on the same footing as the other slow-roll parameters. In order to emphasize the 
difference between generalized kinetic terms and canonical kinetic terms, we have however 
separated the terms proportional ^ from the others. 
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The system of equations (l49ll that we have obtained is of the form 



u" + 2Lu' + Qn = . (53) 

The matrix coefficient for the first order time derivative is 2L = 2E/r, where E is an an- 
tisymmetric matrix, hnear in the slow-roll parameters. Let us introduce a time- dependent 
orthogonal matrix R which satisfies R' = — LR. Note that this is possible only if L is an 
antisymmetric matrix. Reexpressing the above equation ( !53l ) in terms of a new matrix 
vector V, defined hj u = Rf , one can eliminate the terms proportional to the first order 
time derivative and we obtain the following equation 

v" + R"^ (-L^ - L' + Q) Rt; = 0. (54) 

At linear order in the slow-roll parameters, one finds 

- - L' ~ -i^E. (55) 

Therefore, apart from the trivial part proportional to the identity matrix, the combination 
— — L' + Q contains 

-(E + M) - - ^ r^.. - ^sl vss - ^ce{l + 4))^ ^^^^ 
which is a symmetric matrix. 

We now assume that the slow-roll parameters vary sufficiently slowly during the few 
e-folds when the given scale crosses out the Hubble radius. We thus replace the time- 
dependent matrix on the right hand side of (l56l l by the same matrix evaluated at Hubble 
crossing, i.e. for k = aH, and the only remaining time dependence appears in the global 
coefficient 3/t'^. One can now diagonalize this matrix by introducing the time-independent 
rotation matrix 

y smH* cosfc)* J 

so that 

R;1(M + E)R, = Diag(Ai,A2). (58) 

In particular, one can easily compute the following linear combinations, which will be 
useful later: 

Ai + A2 = 3 {r]„a + Vss - 2e - ^ce) , (59) 

(Ai - A2) sin 26, = 6 (r^,, - ^sl) , (60) 

(Ai - A2) cos 26, = 3 [r]^^ - r]ss - 2e + ^ce{l + 2sl)) , (61) 
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where the right hand sides of the three above equations are evaluated at k = aH. 



Similarly, the rotation matrix R is slowly varying per efold, since (dR/d lna)R"^ = E, 
where E is linear in slow-roll parameters. Around Hubble crossing, one can thus replace 
R by its value R* at Hubble crossing. 

By introducing 

w = R;^R,v, (62) 
the system of equations can be written as two independent equations of the form 

1 



(2 + 3Xa] 



WA = 0, (^ = 1,2) 



with 



Defining 



9 

- + 3Xa 



the solution of (l63ll with the proper asymptotic behaviour can be written as 



Wa 



^ i(A*A+l/2)7r/2, 



e-- V-rHj^}i-kT)eAik), 



(63) 
(64) 
(65) 

(66) 



where H^^ is the Hankel function of the first kind of order n and the ca are two indepen- 
dent normalized Gaussian random variables so that 



{cAik)) = 0, {eA{k)eUk')) = SAsS^'^k - k'). 



(67) 



Using the independence of the variables Wi and W2, one can express the correlations 
for the variables u^j and Ug around Hubble crossing time as 

a^{QlQa) = cos^ e^{wlwi) + sin^ e^{wlw2) (68) 

a^iSs^Qa) = ^sm2Q, (^{wlwi) - {wlw2)^ (69) 

0^(53^53) = sm^e{w\wi) + cos^e{wlw2) (70) 
where one can substitute 

(^W) = li-r)\H;^li-kr)\' ^ ^ (^•^^(-/.r). (71) 
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This yields 

Vq^ = 

Vss = 

where we have used 



(^^^ (1 - 2e,) [cos^ Tii-kr) + sin^ r2{-kT 



§1 



sin 20* 



(^^^ (1 - 2e,) [sin^ e, Ti{-kT) + cos^ 6, J'2(-A;r 



1 + e. 



T 



(72) 
(73) 
(74) 

(75) 



At this stage, it is worth noting that our derivation is still valid if the parameter r^^^, which 
corresponds to the curvature of the potential along the direction orthogonal to the field 
trajectory, is not small. In this case, the entropy fiuctuations are effectively suppressed 
and only adiabatic fluctuations are generated. As far as the perturbations are concerned, 
this particular situation is similar to the single field case. 

A further simplification occurs when 1, in which case /xa — | + The functions 

J^a{x) can be expanded as 

J^a{x) = ^x'\H,/,{x)\' (1 + 2XAfix)) = (1 + x^) (1 + 2XAfix)) , (76) 



with 



fix) = Re 



1 dHl!\x) 



Hi%ix) 



(77) 



;i=3/2 , 



Using the relations (15911611) and (1641) . we finally get for the curvature and entropy 
perturbations defined in ([271) and (l28l) . the following expressions 



C-ns 
Vs 



H: 



2 \ 2 



27r6"* 



:i + k'r 



1 - 2e* + (6e* - 2ri„^^ - 2{*s^^C0*) / 
k 



k 



p-\\l + k\') (2^4 - 2r^.,.) / f 
27rcr* / \ 



H 

27rcr* 



2 \ 2 



'l + k'T 



2„2> 



1 - 2e, + (2e, - 2r/,,, + 2^,(1 + SeJcg*) / 



(78) 
(79) 
(80) 



Let us comment these results. First, one can verify that, for canonical kinetic terms 
(i.e. ^ = 0), we recover the results of [22] if we replace the factor (1 + fc^r^) by 1 and 
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the function f^—kr) by the number C = 2 — ln2 — 7 ^ 0.7296, where 7 ^ 0.5772 is 
the Euler-Mascheroni constant. Our final expression depends exphcitly on r and allows 
us a more precise estimate of the spectra around the time of Hubble crossing. As we 
will see explicitly later, there are some inflationary scenarios where the amplitude of the 
curvature perturbation spectrum evolves very quickly after Hubble crossing and never 
reaches its asymptotic value (corresponding to the limit kr 0). In these cases, one 
needs to evaluate more precisely the amplitude around Hubble crossing, which our more 
detailed formula enables to do. 

Another difference with [22] is that we derived the adiabatic and isocurvature spectra 
by working directly with the equations for the adiabatic and entropy components, instead 
of working with the initial scalar flelds. 



4 Evolution of perturbations on super-Hubble scales 

When the isocurvature modes are suppressed, for instance if the effective mass along the 
isocurvature direction is large with respect to the Hubble parameter the flnal adiabatic 
spectrum can be computed simply by taking the usual single-fleld result applied to the 
adiabatic direction: 

where all the quantities are evaluated at Hubble crossing. The simpliflcation works be- 
cause, in this particular situation where isocurvature fluctuations are absent, the curvature 
perturbation remains frozen on super-Hubble scales. 

If isocurvature modes are present however, they will affect the super-Hubble evolution 
of the adiabatic perturbations because they will act as a source term on the right hand 
side of the equation governing the evolution of the curvature perturbation [7j (and [28] 
for the non-linear generalisation). 

In order to obtain the flnal power spectra and correlations, and to compare the pre- 
dictions of a multi-inflaton model with observations, one must then solve the coupled 
system of differential equations (I34tl35p . In general, a numerical approach is necessary 
and will be considered in the next section. In some particular cases, within the slow-roll 
approximation, one can solve analytically the equations of motion on super-Hubble scales. 
We now discuss these cases in the rest of this section. 

Following |21j, we can then write eqs. (l34l) and (l35l) in the slow-roll approximation as: 
~ AHQ^ + BH6s and 6s ~ DH6s , (82) 
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where: 



D 



B 



A 



-T]„a + 2e - ^ceSg 
-riss + ^ce{l + si) 




(83) 
(84) 
(85) 



Qualitatively, it is clear that if the isocurvature perturbations do not decay very fast, there 
is a strong interaction between the adiabatic and isocurvature perturbations, whenever 
the coefficient B becomes large, i.e. when the classical trajectory makes a sharp turn in 
the field space or when ^ is relatively large and infiation is driven at least partially by the 
field X {se 7^ 0). For constant A, B, D, the equations (l82l) can be solved explicitly to give 



Q^{N) ~ e^^Q^, + -(e^^ - e^^)(5s, and 6s{N) ~ e^^5s, (86) 



where N stands for the number of efolds after Hubble crossing. Taking into account that 
{H/&) ~ (//^,/cr*)e~^^, we can express the power spectra and correlations as: 



where 'y = D — A. The quantities V-ji*, Cns* and Vs* correspond to the asymptotic limit, 
i.e. when kr 0, of the expressions fl78tl80p . 

The use of this approximation is in practice rather limited because it relies on the 
assumption that the slow-roll parameters are time-independent between Hubble crossing 
and the final time. In most cases, this approximation, which we call constant slow- 
roll approximation, holds only for a few e-folds and breaks down long before the end of 
infiation. 

In some simple infiationary models, there exists an analytical approach to compute 
analytically the evolution of the perturbations on super-Hubble scales. This is the case 
for double infiation with canonical kinetic terms (i.e. 6 = 0) and potential 



where the equations of motion for the metric perturbation $ and the two scalar field 
perturbations can be integrated explicitly in the slow-roll approximation and on super- 
Hubble scales [5] (this can be seen as a particular case within a more general context 



B 




(90) 
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discussed in [9]). One finds 



$^-^+2C3 \^,^^,,.r, n,.;; . (91) 



~r ~ ^*^3 2 2 I 2^2' ~ — + ^*-'3 2 2 I 2^2' I^^^J 

where Ci and C3 are time-independent constants of integration. 

The curvature and isocurvature perturbations during inflation are respectively 

7^ = $ + /7^^^^±i^, S = H^^2^^4^. (93) 

^2 _^ 02 ^2 _^ 02 



By plugging the solutions f l9Tll92ll into the above expressions, one obtains the explicit 
evolution of the adiabatic and isocurvature perturbations, knowing that the background 
evolution is given by 

(sina)2/(^'-i) 

X = 2Mp^/s sin a, (f) = IMp^l cos a, s = Sq- .^Ba/rBZ n (^4) 

(cos a)"''" 

where s = — ln(a/ae) is the number of e-folds between a given instant and the end of 
inflation, and R = m^/m^i,. 

Note that the isocurvature perturbation S that we have deflned during inflation, fol- 
lowing other works, is proportional but does not coincide with the isocurvature pertur- 
bation S'rad deflned during the radiation era. In the scenario discussed in [10], where the 
heavy held x decays into dark matter, the isocurvature perturbation S'rad = ^cdm — (3/4)5^ 
is related to the perturbations during inflation via the relation 

Another method to calculate the flnal curvature perturbations is the so-called 6N 
formalism [311 [32l [33] . In practice however, this method requires the expression of the 
number of e-folds as a function of the initial values of the scalar flelds and except in a few 
simple cases where this expression can be determined analytically, a numerical approach is 
also needed in the general case. Moreover, the approach we have adopted allows to follow 
not only the evolution of the curvature perturbation but also that of the isocurvature 
perturbation. This is important if some isocurvature perturbations survive after the end 
of inflation. Their evolution then depends on the details of the processes which occur at 
the end inflation and after, in particular the reheating (or preheating) processes, which 
goes beyond the scope of this study. 
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5 Numerical analysis 



In Section [3l we studied the spectra and correlations of the perturbations in the vicinity 
of the Hubble crossing and we obtained analytic approximations (rf8l ) - ([80l) . The aim 
of the present section is to confront these expressions with the result of a numerical 
integration of the equations of motion (i34ll -(l35l). We would also like to study numerically 
the super-Hubble dynamics of the perturbations, which is often the only way to calculate 
the inflationary observables such as the spectral index with a precision required by the 
present and forthcoming observations. 



5.1 Numerical procedure 

Our numerical procedure, similar to that of [23], is the following. In order to take into 
account the statistical independence of the adiabatic and isocurvature perturbations deep 
inside the Hubble radius, we integrate eqs. (l34ll -(l35ll twice: first with the initial value 
of Qa corresponding to the Minkowski-like vacuum and 5s = 0, then with the initial 
value of 5s corresponding to the Minkowski-like vacuum and = 0. Unless stated 
otherwise, we impose the initial conditions 8 efolds before the Hubble crossing. The initial 
conditions also include the slow-roll for the background fields. The evolution proceeds 
along a background trajectory which provides a sufficient number of efolds before the end 
of infiation (50-60, depending on the model). We identify the end of infiation, at which we 
terminate the evolution, when e = 1. As the outcome of the first (second) run, we obtain 
the curvature and entropy perturbations, IZi and Si (7^.2 and S^)- We then calculate the 
spectra and correlations as: 

Vn = ;^(|7^lp + |7^2p) (96) 

Vs = l^d^ir+i^^D (97) 



C 



ns 



n\Si + nls^ . (98) 



27r2 

We shall sometimes describe the correlations using the relative correlation coefficient: 

C = 1 '^'^'^* (99) 

The value of C lies between and 1, and it indicates to what extent the final curvature 
perturbations result from the interactions with the isocurvature perturbations. 



17 



(a) 



-8 (start) 
(Hubble crossing) 



61 (end) 50 40 30 




(c) 



2,6 - 

2.4 - 

2.2 - 

2 - 

1.8 - 

1.6 - 

1.4 - 

1.2 7 

1 r 

0.05 0.1 0.15 0.2 



50 (end) 



-8 (stort) 
(Hubble crossing) 

10 

20 
30 I 
40 / 



0.25 



Figure 1: Examples of classical inflationary trajectories for double inflation with canonical 
kinetic terms (left), double inflation with non-canonical kinetic terms (center) and roulette 
inflation (right). The details of the models are described in Section [57^ Subsequent tens 
of efolds are indicated along the curves. 



5.2 Examples of inflationary models 

There is an enormous number of examples of inflationary models. Here we restrict our 
analysis to just three cases described below. We will use these examples to check the 
analytical results of the preceding Sections and to illustrate some generic features in the 
evolution of adiabatic and isocurvature perturbations. 



5.2.1 Double inflation with canonical kinetic terms 

Double inflation (with 6 = 0) is certainly the most thoroughly studied example of multi- 
fleld inflation. It employs the potential 

V{<j),x) = \my + ^mlx'. (100) 

In order to make definite calculations, we set Tm^ = and we choose the initial condi- 
tions = Xi (later we shall also comment on the case (pi = 50xj), 8 efolds before the scale 
we consider leaves the Hubble radius, after which inflation goes on for about ~ 60 efolds. 
The classical trajectory in field space is shown in Figure [U In our example, the trajectory 
is strongly bent roughly 35th efolds after the Hubble crossing. As we shall see, it is at 
this moment when the adiabatic perturbations are strongly fed by the isocurvature ones. 
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5.2.2 Double inflation with non-canonical kinetic terms 



In order to concentrate just on the effects due to the non-canonical nature of the kinetic 
terms, we consider a very simple generalization of the previous example by taking 6(0) = 
—(p/Mp and m<^ = in (1 1001) . We choose the initial conditions so that (pi = 0. Then it 
is almost exlusively the field x which slides down to the minimum of the potential during 
inflation, but due to non-canonicality of the kinetic terms, the interaction of x with 
drives the latter slightly away from zero. The classical trajectory in the field space is 
shown in Figured! 

5.2.3 Roulette inflation 

Recently, inflation in the large volume compactiflcation scheme in the type IIB string 
theory model has been investigated in [H] (see also f34]). In our notation, this model can 
be effectively described by 



where ip{(p) = {(p/Mp)^/^ and 6o, Vi, A ^ire functions of the parameters of the underlying 
string model. A generic feature of the potential (|102p is that it has an infinite number 
of minima arranged periodically in x ^ind a plateau for large values of (p, admitting a 
large variety of inflationary trajectories, which may end at different minima even if they 
originate from neighboring points in the held space - hence the model has been dubbed 
roulette inflation. In this work, we adopt the parameter set no. 1 (in Planck units: 
bo = -11; Vo = 9.0 X 10"^^ Vi = 3.2 x 10"^ V2 = 1.1 x 10"^; and /?i = 9.4 x 10^ 
/?2 = 27r/3) from [14] and choose the particular inflationary trajectory shown in Figure [H 
For this trajectory, the factor b^Mp is rather large, of the order 10^, but the effect of the 
non-canonical kinetic terms is strongly suppressed by a very small value of e on the plateau 
of the potential. The smallness of e also suppresses the energy scale of inflation and one 
needs a smaller number of efolds than in the models described above. For deflniteness, we 
assumed that there are ~ 50 efolds between the moment that the scale of interest crosses 
the Hubble radius and the end of inflation. 

5.3 Numerical results for the perturbations 



For the three inflationary models described in Section 15.21 we performed the numerical 
analysis, as described in Section 15. 1[ Here, we discuss the outcome for the spectra and 




(101) 



and 



V{(P, x) = Vo + V.^/ii^e-'"^^^'^^ + V2i:{(P)e-'''^'^^'> cos(/?2X) 



(102) 
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correlations in Figures [SHU In the right panel of each Figure we plot the evolution of Vn 
and Vs^ normalized to the single-field result given in eq. (ISTll . as well as the evolution of 
the correlation coefficient C defined in (l99l l and the parameter -B, defined in (l84l) . which is 
the coupling between the isocurvature and the curvature perturbations. These quantities 
are plotted as functions of the number of efolds N after Hubble crossing. Left panels of 
each Figures [2H4] are basically close-ups of the right ones to the vicinity of the Hubble 
crossing. There, we plot the evolution of Vn, "Ps and C-jiSi normalized to the single- 
field result given in eq. flSTj) . These are shown as functions of the variable {k/aH)~^ 
which allows us to compare directly the numerical results with the predictions of the eqs. 
(I78ll -(l80l). Note that in the leading order in the slow-roll parameters ln{aH/k) = N, 
hence the logarithmic scale in the left panels directly corresponds to the linear scale in 
the right panels. In Figures [IHH we also plot the evolution of the spectra, denoted by 
the superscript (a), when one assumes the constant slow-roll approximation after Hubble 
crossing, i.e. when one uses eqs. (I87l ). ( l88l) and f l89l) . 

For completeness, we also discuss briefiy the particular case where the generation of 
isocurvature modes is effectively suppressed, situation which applies to some of the models 
discussed in the literature for specific parameters and/or initial conditions. 

5.3.1 Double inflation with canonical kinetic terms 

In this example, the field x initially dominates the energy density of the Universe and 
drives the first part of infiation, and only when it is almost at its minimum, infiation is 
further driven by 0. All the slow-roll parameters are small at the Hubble crossing, which 
makes eqs. (TfHl i - IISOl) an excellent approximation of the numerical solutions of the equations 
of motion. Due to the smallness of -B = —2rjcjs ~ 2d6'/dA^ right after the Hubble crossing, 
the curvature perturbations become practically constant during the x-domination. At the 
transition to 0-dominated inflation B becomes large, which leads to a sizable increment 
in Vn because of interaction with the isocurvature perturbations. During 0-dominated 
inflation the trajectory is almost straight again, the isocurvature perturbations decay 
quickly and the curvature perturbations are frozen at the value acquired at the transition. 

This model has the advantage that the numerical results can be directly compared to 
analytical ones, as the time evolution of the perturbations can be solved without assuming 
constancy of the slow-roll parameters [10], and we find a good agreement between two 
approaches. 
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Figure 2: Predictions for the spectra and correlations of the perturbations in double in- 
flation with canonical kinetic terms. Thick lines show the numerical results for V-ji, Vs 
and C-jis or C normalized to the single-field result ( fgljj . respectively. Circles, stars and 
squares indicate the predictions of eqs. ( TT^) . ( TT^J and (E^, respectively. Thin dashed lines 
indicate the predictions of eqs. ( fg?IJ . ( fggjJ and ( fgPIJ . respectively. The coupling B between 
the curvature and isocurvature perturbations is also shown. 



5.3.2 Double inflation with non-canonical kinetic terms 

In this example, the background trajectory is almost straight. However, the slow-roll 
parameter e is around 0.1, which makes the coupling B large throughout the entire infla- 
tionary era. Figure [3] shows that eqs. (i78l) -(l80l) are a good approximation for the spectra 
and correlations at the Hubble crossing, k/aH = 1, but it is no longer true at super- 
Hubble scales, k/aH = 0.1 or 0.01, because the isocurvature perturbations already start 
feeding the curvature ones sizably. As a result, the flnal curvature perturbations originate 
almost exclusively from the interactions with the isocurvature modes, not from the fluc- 
tuations along the inflationary trajectory, which makes the relative correlation coefficient 
very close to 1. 
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Figure 3: Predictions for the spectra and correlations of the perturbations in double in- 
flation with non- canonical kinetic terms. Thick lines show the numerical results for V-ji, 
Vs and C-jis or C normalized to the single-field result ( fgljj . respectively. Circles, stars and 
squares indicate the predictions of eqs. ( TT^) . ( TTglJ and (E^, respectively. Thin dashed lines 
indicate the predictions of eqs. ( fg?IJ . ( fggjJ and ( fgPIJ . respectively. The coupling B between 
the curvature and isocurvature perturbations is also shown. 



5.3.3 Roulette inflation 

As we already argued in Section [521 most of the inflationary trajectory in this example lies 
on the plateau of the potential (I102p . the slow-roll parameter e is very small, which makes 
the direct impact of the non-canonicality negligible. The trajectory is, however, strongly 
curved in the field space and the interaction between the isocurvature and curvature modes 
is still important. Again, eqs. (rf8l) - ([80l) accurately predict the spectra and correlations in 
the vicinity of the Hubble crossing, with deviations on super-Hubble scales resulting from 
the sourcing of the curvature perturbations by the isocurvature ones. Eventually, most of 
the curvature perturbations arise through this effect. 
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Figure 4: Predictions for the spectra and correlations of the perturbations in roulette 
inflation. Thick lines show the numerical results for V-ji, o-nd C-jis or C normalized to 
the single-field result ( fgljj . respectively. Circles, stars and squares indicate the predictions 
of eqs. fM) . ( TTPI) and ^ES), respectively. Thin dashed lines indicate the predictions of eqs. 
(E3j, ^SR) and ^EM), respectively. The coupling B between the curvature and isocurvature 
perturbations is also shown. 



5.3.4 Impact of the non-canonical terms 

In order to estimate the impact of the non canonical kinetic terms, one can separate each of 
the coefficients fl36tl39p into two parts: the terms that depend explicitly on the derivatives 
of non-canonical function h and those which do not. For example, the coefficient that is 
directly responsible for the transfer of isocurvature modes into curvature modes, C^s is 
decomposed into a 'canonical' component C'^} and a 'non-canonical' component C^f* : 

c^, = &:} + c^:f\ (103) 

with 

C^^:^ = 6i/f + ^ + 2K. + ^, Ct'=2h,{slV^-clV,). (104) 
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Figure 5: Evolution of the coefficients Ca] and Cj^^^ defined in Iil04\ ), parametrizing the 
coupling between the curvature and isocurvature modes: (a) for double inflation with non- 
canonical kinetic terms and (b) for roulette inflation. 

The other coefficients can be decomposed similarljQ. 

For the two models with non-canonical kinetic terms, we have compared in Fig. [5] the 
canonical and non-canonical contributions of the various coefficients. For double inflation 
with non-canonical kinetic terms, the non-canonical contribution in C^s is dominant and 
therefore plays a crucial role in the evolution of the curvature mode. For roulette inflation, 
as already noticed earlier, the non-canonical contributions turn out to be negligible. 

5.3.5 Effectively single-field cases 

Many supergravity- or string-inspired models aim at describing supersymmetry breaking 
and inflation in a unified framework. Often, despite the presence of many scalar fields, one 
can find model parameters and initial conditions such that only for one combination of the 

^Note that there is some arbitrariness in this decomposition: the decomposition would be different if 
one expresses Va in terms of 9 via the background equation (|33l) . 
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Figure 6: Predictions for the spectra and correlations of the perturbations in the effective 
single-field case. The lines show the numerical results for V-jz and Vs, circles correspond 
to prediction of eq. [7^) and stars correspond to the analytic calculation outlined in Section 

no:i 



fields a potential is sufficiently flat to support inflation, e.g. in pseudo-Goldstone inflation 
[TS] . "better racetrack" scenario [Hj, no-scale supergravity models with moduli stabilised 
through D-terms [17] or in D-term uplifted supergravity of In these works, small 

values of the field velocity (i.e. e -C 1) have been ensured by setting the initial conditions 
in the vicinity of the saddle point of the potential, while small curvature of the potential 
along one direction has been obtained by a fine-tuning of the parameters. It has been 
assumed that the isocurvature perturbations decay fast after the Hubble radius crossing 
and do not affect the curvature perturbations even though the trajectory in the field space 
can have sharp turns at later stages of the inflationary evolution. Consequently, the single 
field approximation (ISTl) has been used in these works to account for the spectra of the 
curvature perturbation. 

In the two-field language developed in the present paper, the situation described above 
corresponds to rjss > 0{1) ^ IVeral, |^o-s|,e- This justifies setting 0* = in eqs. (l68l) - (j70l) . 
which is equivalent to the assumption that the curvature and isocurvature modes evolve 
independently. While we can still apply the slow-roll expansion that led to eq. (l78l l for the 



25 



power spectrum of the curvature perturbations, we now have to use the full result (l66ll to 
describe the spectrum of the isocurvature modes. For rjss ^ C'(l), the latter result decays 
very fast for A;|r| ^0 and we conclude that the isocurvature modes become irrelevant for 
the evolution of the curvature perturbations soon after the Hubble radius crossing. 

We checked numerically that the above conclusion applies for the models described 
in [17], which are easily reduced to the two-field case and their inflationary trajectories 
are curved in the field space. For simplicity, we would like, however, to illustrate this 
point with the model of double inflation with standard kinetic terms, described in Section 
15.2.11 for which we set the initial condition (pi = 50xi- Then the heavy held x contributes 
negligibly to the potential energy and rjss* — 0.4. In Figure [H we plot the numerically 
calculated spectra of the curvature and isocurvature perturbations and compare them with 
the analytic approximations outlined above. The decay of the isocurvature modes agrees 
with the solution (l66ll . for which we show two cases: the small solid stars correspond to 
the constant value oirjss, whereas the large empty stars show the result corresponding to 
adjusting the index of the Hankel function in eq. (l66ll to the value of r]ss at a given instant, 
i.e. r]ss = 0.40, 0.42, 0.44 for {k/aH)~^ = 1, 10, 100, respectively. With the isocurvature 
modes absent, the curvature perturbations are excellently described by the single-fleld 
result, which justifles the use of the single-fleld approximations in the situations described 
above. 

5.3.6 Closing discussion 

The three examples presented here show that in multi-fleld inflationary models, a large 
part of the curvature perturbations can originate from interactions between the curvature 
and isocurvature perturbations on super-Hubble scales, not only from quantum fluctua- 
tions along the trajectory at the Hubble exit. In such cases, the single-fleld result (ISTll 
does not provide a correct prediction either for the normalization of the power spectrum 
or for its spectral index 

ns = l + dlnVn/dlnk. (105) 

There are techniques which allow relating the spectral index of the curvature perturba- 
tions. Us, to the spectral indices of the entropy perturbations and the curvature-entropy 
correlations through a set of consistency relations [El [22l [T3j, but all these quantities 
separately depend on the super-Hubble evoulution of the perturbations. Again, we flnd it 
the most straightforward to calculate the spectral index Ug for each model numerically. In 
Table [T], we compare naive estimate Ug ~ l — 6e^ + 2ri„^^ and the predictions of eq. flSTj) for 
the spectral index Ug with the numerical results of Section [531 In our three examples, one 
can see that the single-fleld result signiflcantly overestimates the correct spectral index. 
The discrepancies that we flnd follow from the fact that the two types of perturbations 
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1 — 6e* + 2?7crcr* single-field result full result 


double inflation (canonical) 
double inflation (non-canonical) 
roulette inflation 


0.929 0.982 0.967 
0.953 0.968 0.934 
1.017 1.019 0.932 



Table 1: A comparison between the predictions for the spectral index Ug in the three 
examples of inflationary models described in Section lSTR The third column contains result 
derived from the single-field approximation ( fgljj ; the result of full numerical calculations 
are shown in the fourth column. 

experience the slow-roll of the background fields and the curvature of the infiationary 
potential in a different way. Then, if the final curvature perturbations originate mainly 
from the isocurvature ones, they inherit the features of the isocurvature power spectra at 
the Hubble crossing. 

6 Conclusion 

In this paper, we have studied two-field inflation and extended several previous results on 
curvature and isocurvature perturbations to the case of non-standard kinetic terms. 

First, we have calculated analytically the curvature and isocurvature spectra, as well as 
the correlation, just after Hubble crossing for two-field infiation models, including next-to- 
leading order corrections in the slow-roll approximation. Our results fl78l)-fl80ll generalize 
those of Byrnes and Wands, who assumed only standard kinetic terms. We have also given 
a refined analytical treatment of the spectra around Hubble crossing, which is important 
when the perturbations still evolve after Hubble crossing. 

Second, we have studied numerically the evolution of the curvature and isocurvature 
perturbations after the Hubble crossing. This type of analysis is important since in multi- 
field infiation, in contrast with single-field infiation, the curvature perturbation spectrum 
after infiation is in general different from the curvature perturbation spectra at Hubble 
crossing because of isocurvature perturbations, as first emphasized in |6]. This well- 
known result applies to multi-field infiation with either standard or non-standard kinetic 
terms. This effect has been studied numerically for standard kinetic terms, using the 
decomposition into instantaneous curvature and isocurvature, in the analysis of |23| . We 
have done a similar analysis here for non-standard kinetic terms. In particular, we have 
compared the numerical evolution of the perturbations with an analytical approximation, 
which we denoted the constant slow-roll approximation: this approximation assumes not 
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only that the slow-roll approximation is valid, but also that the slow-roll parameters 
remain almost constant during the subsequent evolution where isocurvature perturbations 
are significant. This approximation has been used in [2T] to compute the "final" spectra in 
the context of non-standard kinetic terms. In most cases, this approximation is however 
not very realistic as we clearly show in our numerical study. We have also estimated 
the impact of non-standard kinetic terms on the coupled evolution of the curvature and 
isocurvature modes. 

There has been a recent interest in constructing inflationary models in the context 
of string theory. These models naturally lead to scalar fields with non-standard kinetic 
terms, to which our analysis can apply. In this work, we have studied a very recent model, 
called 'roulette" inflation, and computed quantitatively the final curvature spectrum, thus 
showing that the infiuence of isocurvature modes on super-Hubble scales turns out to be 
very important. Note that the authors of [Hj computed the final curvature spectrum 
by using a "single inflaton approximation". They stress however that isocurvature effects 
"could produce big effects", which we indeed confirm in our analysis. 

As a message of caution for the readers who are not familiar with the effect of isocurva- 
ture modes in multi-field infiation, we have also computed the spectral index for curvature 
perturbations in a few models and contrasted it with the value that one would naively 
obtain by using the curvature perturbation at Hubble radius like in single field infiation. 

Finally, an interesting question, which goes beyond the scope of this paper, would be 
to investigate in which circumstances these multi-field models could produce isocurvature 
perturbations, after infiation and the reheating phase. These "primordial" isocurvature 
perturbations are today severely constrained by CMB data. 
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